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Damping  Characteristics  of  Fiber  Composites*'*^ 

by 

Zvi  Hashin*^**^ 

Abstract 


Analytical  results  for  complex  moduli  of  uniaxially  fiber  reinforced  mat¬ 
erials  made  of  viscoelastic  matrix  and  elastic  fibers  are  reviewed.  A  general 
method  is  established  to  predict  complex  moduli  and  loss  tangents  of  viscoelas¬ 
tic  laminates  made  of  uniaxially  reinforced  laminae.  Application  of  results 
is  demonstrated  by  two  examples  of  analysis  of  damping  of  structural  vibrations: 
Attenuation  of  vibrations  in  uniaxially  reinforced  Timoshenko  beam  and  torsional 
vibrations  of  laminated  cylinder. 


f*) 

Presented  at  8th  International  Council  of  the  Aeronautical  Sciences  (ICAS) 
Congress,  Amsterdam,  August  28  -  September  2,  1972. 

(**) 


Professor,  Department  of  Materials  Engineering, 
Technion  -  Israel  Institute  of  Technology. 
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1 .  Introduction 

The  ever  increasing  use  of  fiber  composites  for  aero/space  structures 
requires  the  development  of  rational  methods  for  prediction  of  their  rele¬ 
vant  properties  within  engineering  accuracy. 

In  most  current  fiber  composites  the  matrix  is  a  polymer  such  as  epoxy. 
It  is  well  known  that  such  polymers  exhibit  the  effect  of  vibration  damping. 
Therefore  such  damping  effects  will  also  occur  in  composites  in  which  the 
matrix  is  polymeric.  Since  aero/space  structures  are  subjected  to  severe 
vibrational  environment  and  since  vibration  damping  is  beneficial,  the 
quantitative  prediction  of  such  damping  is  of  considerable  engineering  imp¬ 
ortance. 

It  is  of  interest  to  emphasize  the  unique  combination  of  desirable  pro¬ 
perties  which  are  exhibited  by  fiber  composites:  Superior  strength  and 
stiffness,  low  weight  and  vibration  damping.  No  other  materials  seem  to 
possess  this  many  advantages. 

In  order  to  handle  the  problem  analytically,  it  is  assumed  that  the 
matrix  is  linearly  viscoelastic.  Its  dynamic  viscoelastic  properties  can 
then  be  characterized  in  terms  of  the  usual  complex  moduli  of  viscoelasticity 
which  are  assumed  to  be  known  on  the  basis  of  experiments.  The  fibers  are 
represented  as  linear  elastic.  The  composite  with  such  constituents  behaves 
macroscopically  as  a  linear  viscoelastic  body  which  is  characterized  by 
effective  complex  moduli.  There  arise  three  classes  of  important  investi¬ 
gation: 
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(a)  Prediction  of  effective  complex  moduli  of  a  uniaxially  reinforced 
material  on  the  basis  of  matrix  complex  moduli;  fiber  elastic  moduli  and  in¬ 
ternal  geometrical  parameters  such  as  constituent  volume  fractions,  fiber 
shapes,  etc. 

(b)  Prediction  of  the  effective  complex  moduli  of  a  laminate,  whose 
laminae  are  composed  of  uniaxially  reinforced  material,  on  the  basis  of  the 
uniaxial  material  effective  complex  moduli  found  in  (a)  and  the  laminate 
internal  geometry. 

(c)  Viscoelastic  vibration  analysis  of  structures  made  of  fiber  compo¬ 
sites. 

These  different  kinds  of  problems  will  be  discussed  consecutively. 
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Complex  Moduli  of  Uniaxially  Fiber  Reinforced  Materials 

A  general  theory  of  prediction  of  effective  complex  moduli  of  composites 
with  linear  viscoelastic  constituents  has  been  given  previously  [1] ,  [2] ,  [3] . 
It  will  here  suffice  to  discuss  without  proof  some  results  which  are  pertinent 
for  the  present  investigation. 

Let  the  local  average  strains  and  stresses  in  a  composite  be  of  oscillatory 
nature .  Thus : 


a. . 


(2.1) 


where  overbars  denote  average,  x  =  to  is  frequency,  t  is  time  and  latin 
subscripts  range  over  1,  2,  3.  The  effective  complex  moduli  of  a  generally 

anisotropic  composite  are  defined  by  the  relation: 


(2.2) 


where  superscripts  R  and  I  denote  real  and  imaginary  parts  respectively. 

The  assumption  is  made  that  the  fiber  reinforced  material  under  considera¬ 
tion  is  macroscopically  transversely  isotropic  with  respect  to  fiber  direction. 
Then  (2.2a)  assumes  the  form; 
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5il  =  Cii  Ell  +  Ci2  £22  +  Ci2  £33 


^22  =  Ci2  Ell  +  C22  £22  +  C22  £ 


3  3 


033  =  C12  Ell  +  C23  £22  +  C22  £33 


O12  -  2C44  £12 


023  =  (C22  -  C23)£23 


031  =  2C44  £31 


(2.3) 


where  xi  is  in  fiber  direction  and  X2,  X3  are  in  the  transverse  plane.  Fig.  1. 
In  another  notation,  the  complex  moduli  in  (2.3)  are  written: 

Cu  =  n 

C12  =  ^ 

C22  =  k  +  Gj  (2.4) 

C23  =  k  -  G.J, 


Here  k  is  a  transverse  complex  bulk  modulus,  G.p  -  transverse  complex  shear 
modulus  in  X2,  X3  plane  and  G^  -  axial  complex  shear  modulus  in  planes  contain¬ 
ing  fiber  direction  Xi.  The  physical  interpretation  of  n  and  £  is  here  of 
little  interest. 

Inversion  of  (2.3)  is  written  in  the  form: 


1  ~  ~  \  ~ 
En  =  TT— <7ii  -  r— O22  -  r— 033 

EA  EA 


^A  ~  1  ~  ~ 

£22=-  +  5— <722  -  5—033 

Ea  Et  Et 


£33--  ^011  -  ;r-  022  •*•  ^=—033 
E^  Ex  Ex 


(2.5) 


where  E^  and  E.^.  are  complex  Youngs'  moduli  in  axial  (fiber)  direction  and 
transverse  (to  fiber)  direction,  respectively,  and  and  v.p  are  associated 
complex  Poisson's  ratios. 

Establishment  of  analytical  expressions  for  the  various  effective 
complex  moduli  listed  above  in  terms  of  matrix  complex  moduli,  fiber  elastic 
moduli  and  phase  geometry  is  based  on  a  correspondence  principle  [1],  [2] 
which  states:  The  effective  complex  moduli  of  a  viscoelastic  composite  are 
obtained  by  replacement  of  phase  elastic  moduli  by  phase  complex  moduli 
of  a  composite  with  identical  phase  geometry. 

In  the  usual  fiber  reinforced  materials  the  following  conditions  are 
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usually  fulfilled  with  sufficient  accuracy 

(a)  Fibers  are  by  an  order  of  magnitude  stiffer  than  matrix 

(b)  The  matrix  is  isotropic  and  is  viscoelastic  in  shear  only 

(c)  The  matrix  shear  loss  tangent  is  not  larger  than  0.1„  Thus; 


.  (.ill 

tan  = 


<  0.1 
m 


(2.6) 


Under  these  conditions  the  following  results  have  been  shown  [3] 
to  be  valid. 

(d)  The  imaginary  parts  of  the  effective  complex  moduli,  n,  1,  ic, 
niuch  smaller  than  the  imaginary  parts  of  the  effective  shear 

moduli  G.^.,  and  of  E.j,. 

(e)  To  obtain  real  parts  of  all  effective  complex  moduli  it  is 
merely  necessary  to  take  corresponding  expressions  for  effective  elastic 
elastic  moduli  and  to  replace  in  them  matrix  elastic  moduli  by  real 
parts  of  matrix  complex  moduli. 

Some  simple  general  results  which  are  valid  under  the  conditions 
listed  above  will  now  be  given: 

=  \(>U)  \  *  Ej  Vj 


E^(a))  =  ♦  Ej  Vj 

e'c«.)  -  e'm 


(2.7) 


tan  6p  =  =  - — — — —  «  tan  5^ 

"a 
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Here  m  and  f  denote  matrix  and  fibers,  respectively,  and  v  stands  for 
volume  fraction,  tan  5^  is  the  loss  tangent  for  uniaxial  stressing  in  fiber 
direction  while  tan  6^  is  the  corresponding  loss  tangent  for  the  isotropic 
matrix . 


For  axial  shear 


1  +  V, 


,R  ^  *  ''f 


G^Cto)  =  G*'  ^ 
m  1  - 


I  1  ^  Vf 


GtCw)  =  G^  T - 

A"-  •'  m  1  -  V 


(2.8) 


For  transverse  shear 


tan  6  =  tan  6’!l! 

Ga  G 


^  x.m 


tan  6„  =  tan  6^ 

U 


(2.9) 


An  expression  for  G^  has  been  given  in  [2] ,  [3] .  Results  for  other 
complex  moduli  may  also  be  found  in  these  references. 
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3.  Complex  Moduli  of  Laminates 


The  laminates  to  be  considered  are  composed  of  plane  laminae  of 
uniaxially  fiber  reinforced  materials,  the  direction  of  reinforcement 
being  different  in  each  lamina.  The  laminate  is  referred  to  a  fixed 
coordinate  system  xi,  xz,  xt  where  xi,X2  are  in  the  plane  of  the  lami¬ 
nae  and  X3  is  normal  to  it.  Fig.  2.  The  nth  lamina  in  the  laminate  is 
referred  to  a  material  system  of  axes  X3  where  xi^^^  is 


normal  to  the  fibers  and  X3  coincides  with  the  laminate  X3.  The  posi- 
Cnl  fnl 

tion  of  the  xj  ,  xz  system  is  defined  with  respect  to  the  Xi,  xa 
system  by  the  reinforcement  angle 


Qjj  =  3c  (xi^"\  Xi) 


(3.1) 


Fundamental  assumptions  of  fiber  composite  laminate  theory  are:  (a)  Any 
lamina  can  be  replaced  by  a  homogenous  material  whose  properties  are  the 
effective  properties  of  the  uniaxial  FRM  of  which  the  lamina  is  made. 

(b)  The  laminae  are  in  statesof  plane  stress. 

First  an  elastic  laminate  will  be  considered.  The  plane  stress-strain 
relations  of  a  lamina  referred  to  its  material  system  of  axes  xi^^\  X2 


are  then; 


^  (n)  1  „  (n)  ''A  „  (n) 

Ell  -^022 


1  rr  (n) 
g-  022" 


(3.2) 
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where : 

-  axial  Youngs  modulus  (in  fiber  direction) 

-  associated  axial  Poisson’s  ratio 

-  transverse  Youngs  modulus  (normal  to  fibers) 

-  axial  shear  modulus  (in  x/"^,  X2^"^  plane) 

The  inverse  of  (3.2)  is: 


=  Cii  +  Ci2  £22^"^ 


a22^’^^  =  Ci2  £11^”^  +  C22  £22^”^ 


(3.3) 


where 


012^^^  =  2Ci»i4  £12 


Cl  1  =  Cl 


Ci2  -  Ci , 


1  - 

E.J,  A 


Vt 


(3.4) 


C22  =  C2222  = 


A 

1“  ^ 
T 


Citi*  =  C1212  =  G. 
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In  terms  of  the  four  index  moduli  in  (3.4),  (3.3)  can  be  written 


compactly  as: 


(3.5) 


where  here  and  from  now  on  Greek  indices  range  over  1,2. 

For  the  sake  of  simplicity  there  will  be  considered  the  special 
group  of  laminates  in  which  the  application  of  membrane  force  in  the 
plane  of  the  laminate  does  not  induce  bending  or  torsion  in  the  laminae. 
The  most  important  kind  of  laminate  which  fulfills  this  requirement  is  a 
symmetric  laminate.  Such  a  laminate  has  the  property  that  its  middle 
plane  is  a  plane  of  symmetry  for  the  geometry  and  elastic  moduli  of  the 
laminate.  The  laminate  is  thus  composed  of  laminae  pairs  in  each  of 
which  the  laminae  are  of  same  thickness,  are  symmetrically  located  with 
respect  to  the  middle  plane  and  have  the  same  elastic  properties  with 
respect  to  the  xi,  X2  system. 

The  last  condition  is  most  commonly  fulfilled  by  laminae  made  of 
identical  material  and  same  reinforcement  angle  0^^  (3.1),  in  each  pair. 
The  elastic  stress-strain  relation  of  such  a  laminate  is  given  by 


No 

ap  p*  -  r„>i 


(3.6) 


where 


-  applied  average  plane  stress 

-  average  strain 

h  -  laminate  thickness 

elastic  moduli  of  laminate 

tj^  -  thickness  of  n^  lamina 

N  -  number  of  laminae 

^aByS  ~  laminae  elastic  moduli  (3.4)  transformed  to  the  xi,  X2 
system  of  axes. 


For  establishment  of  the  result  (3.6)  see  e,g.  [4].  Proof  that  (3.6) 
is  based  on  an  elasticity  solution  which  is  exact  in  the  Saint  Venant  sense 
for  a  laminate  whose  thickness  is  small  compared  to  its  plane  dimensions 
has  been  given  by  B.W.  Rosen  (unpublished) . 

Let  it  be  now  assimed  that  the  laminate  is  viscoelastic  but  remains 


symmetric  as  described  above.  The  laminate  is  subjected  to  oscillatory 
membrane  loads 


^3  "  ^aB  ® 


(3.7) 


The  average  stresses  associated  with  (3.7)  are  then: 


N  N 

-  aB  aB  ^iwt  =  „ia)t 
‘’oS  -  -K - IT  ®  =  “aB  ® 


(3.8) 


The  strain  response  of  the  laminate  is 
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(3.9) 


The  relation  between  and  is  written; 


(3.10) 


where  C*^  g  are  the  effective  complex  moduli  of  the  laminate, 
otpyo 

It  follows  by  the  general  correspondence  principle  of  [1]  which  was 
quoted  above  that  can  be  expressed  in  form  (3.6b).  Thus: 


n=i 


(3.11) 


The  single  laminae  complex  moduli  in  (3.11)  are  interpreted  as 

follows:  In  the  material  axes  xa^^^  of  the  nth  lamina  the  stresses 

and  strains  are: 


a 

~  (n)  lut 

=  ® 

(3.12) 

e 

^  g  (n)  ^ia)t 

aB 

The  relation  between 

j  ~  (n) 

in  (3.12)  is  of  type 

(3. 3-3. 4) . 

Thus: 

^  14 
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By  tensor  transformation: 

^  ^Ciiu  =  =  Ciicos**0  +  C22sin‘*0  + 

n  n 

+  2Ci2Cos^0  sin*0^  +  4C44Cos^0^sin^0 
n  n  tt  n 

^  ^Cii22  =  ^^^Ci2  =  (Cu  +  C22)cos^0  sin^0  + 

n  n 

+  Ci2(cos‘*0  +  sin‘*0^)  -  4C44Cos^0  sin^0 
n  n  n  n 

^  ^C2222  =  ^*^^022  =  Ciisin'*0  +  C22COS*^0  + 

n  n 

+  2Ci2Cos^0  sin^0  +  4C44COS^0  sin^0 
n  ’  n  n 

(3.16) 

^  ^Cui2  =  ^^^Ci4  =  -  Ciicos®0  sin  0  + 

n  n 

+  C22COS0  sin^6  +  Ci2(cos^0  sin0„  -  cos0  sin®0  )  + 
n  n  n  n  n  n*^ 

+  2C44 (cos^0^sin0  -  COS0  sin®0  ) 
n  n  n  n 

^^^£2212=  ^”^£24  =  -  Cucos0  sin^0  + 

n  n 

+  C22cos^0  sin0  +  Ci2(cos0  sin®0  -  cos®0  sin0  )  + 
n  n  n  n  n  n 

+  2C44(cos0  sin^0_  -  cos^0  sin0n) 
n  n  n  " 
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^"^Ci212  =  =  (Cll  +  C22)cos2e^sin20^  - 

-  2Ci2Cos^e^sin^9^  +  C4*,(cos^0^  -  sin^6^)^ 


(3ol6 
cont ’ d) 


The  preceding  developments  together  with  the  results  for  complex  modu¬ 
li  of  uniaxial ly  fiber  reinforced  materials  define  the  computation  methods 
of  the  effective  complex  moduli  of  symmetric  laminates  as  expressed  by  (3,11) 
For  practical  purposes  it  is  frequently  necessary  to  compute  the  effec¬ 
tive  complex  compliances  which  are  defined  as  the  inverse  of  (3,11),  In 
matrix  notation: 

C*  •  S*  =  J  (3.17) 

where  S*  denotes  the  effective  complex  compliance  matrix  and  J  is  the  unit 
matrix.  Separation  of  (3.17)  into  real  and  imaginary  parts  yields 

,  S*^  =  J  (a) 

(3,18) 

R  T  T  R 

c*  •  S*  C*  »  S*  =  0 


Great  facilitation  is  achieved  if  it  is  noted  that  in  view  of  (2.6)  the 
second  term  in  the  left  side  of  (3.18a).  can  be  neglected  with  respect  to 
the  first.  It  then  follows: 


S* 


R=  j 


(a) 

(3,19) 

(b) 


R  I 

Thus,  once  C*  and  £*  have  been  computed  from  (3.11),  (3.19)  define  the 

effective  complex  compliance  matrix  by  simple  real  matrix  operations. 

Another  important  simplification  is  obtained  if  in  the  laminate 

to  each  pair  with  reinforcement  angle  6^  and  thickness  t^  corresponds 

another  pair  with  reinforcement  -6^  and  same  thickness  t  .  It  is  then 

n  n 

easily  realized  by  the  form  of  (3.16)  that  all  contributions  to  (3.11) 
of  terms  with  odd  powers  of  cosG^  and  sinS^  cancel  mutually.  Thus  in 
this  event  the  effective  complex  moduli  matrix  (3.11)  has  the  form 


Cnii  Ch22  0 

=  Cil22  C2222  0 

0  0  C1212 


(3.20) 


and  so  the  laminate  is  macroscopically  orthotropic.  The  situation  just 

described  is  of  frequent  practical  occurrence.  For  example  a  symmetric 

laminate  with  laminae  reinforcement  in  0  =  0,90®,  ±  45®  directions. 

n 
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4.  Structural  Applications 

4.1  Free  flexural  vibrations  of  a 
fiber  reinforced  beam 


As  a  first  example  there  is  considered  the  case  of  free  flexural  vibrations 
of  a  simply  supported  beam  which  is  uniaxial ly  reinforced  in  beam  axis  direction. 
The  purpose  of  the  investigation  is  to  compare  vibration  damping  due  to  matrix 
viscoelasticity  on  the  basis  of  the  usual  theory  which  neglects  the  effect  of 
shear  and  on  the  basis  of  the  more  refined  Timoshenko  theory  which  takes  into  ac¬ 
count  shear  as  well  as  rotatory  inertia.  For  isotropic  materials,  in  which  the 
complex  Young's  modulus  loss  tangent  and  the  complex  shear  modulus  loss  tangent 
are  of  same  order,  the  added  effect  of  shear  and  rotatory  inertia  is  small  for 
vibration  modes  of  low  order  and  for  long  beams.  In  the  present  case,  however, 
where  the  axial  Young's  modulus  loss  tangent  is  by  an  order  of  magnitude  smaller 
than  that  of  the  axial  shear  modulus  (section  2)  the  situation  is  quite  different 
as  will  be  shown  below: 

Considering  only  the  effect  of  flexure  the  differential  equation  of  the 
freely  vibrating  beam  is: 


2  3  w 

C  -K— rr  + 


3^w 


Sw"  3t2 


=  0 


c‘  = 


pA 


(a) 

(b) 


(4.1) 
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where : 

-  complex  axial  Young's  modulus 

I  -  moment  of  inertia 
A  -  area  of  cross  section 
p  -  density 

Boundary  conditions  of  free  support  are: 

3^w 

W,  ^  =0  X  =  o,Jl  (4.2) 

Conventional  viscoelastic  vibrations  analysis  shows  that  the  modes 
of  vibration  are  given  by 

WjjCx,t)  =  \sin  exp(-j  w^tan6gt)  exp(iw^t) 

C4.3) 

(4.4) 


Equ.  (4.4b)  is  a  repetition  of  (2.7d).  The  results  are  valid  for  small 

enough  loss  tangents,  of  order  (2.6).  The  attenuation  n  is  defined  by 

n  ' 

/ 

\  ^  tandg 


(4.5) 
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Next  the  same  beam  is  considered  in  Timoshenko  fashion,  with  shear 
and  rotatory  inertia.  By  the  correspondence  principle  for  viscoelastic 
vibrations  [5] ,  the  elastic  Timoshenko  beam  equation  [6]  transforris  into 


„2  3‘*w  ^  .  kr^p  9‘*w  _ 


(4.6) 


where  c^  is  given  by  (4.1b),  is  the  axial  complex  shear  modulus  (2.8a), 
k  is  the  strength  of  materials  shear  shape  factor  of  the  section  and 
r^  =  I/A. 

Equ.  (4.6)  with  boundary  conditions  of  simple  support  admits  a  solution 
of  the  form,  [6] : 

w  (x,t)  =  A  sin  —  e^V  (4.7) 

n  n  £ 

where  is  a  solution  of  the  complex  frequency  equation: 

kE.  ,2 

c^a  [1+a  2r2(i  +  =  0  (a) 

"  ”  Ga  "  Ga  " 


(4.9) 


It  is  seen  that  (4.10)  is  the  frequency  equation  of  an  elastic 
Timoshenko  beam  in  terms  of  real  parts  of  complex  moduli.  Its  validity  is 
based  on  the  usual  additional  assumption  that  real  parts  of  complex  moduli 
vary  sufficiently  slowly  with  frequency.  Once  has  been  computed  from 
(4.10),  can  be  computed  from  (4.11). 

A  mostly  sufficient  accurate  approximation  for  is: 


For  slender  beams  and  low  modes  the  first  term  in  each  of  numerator 
and  denominator  of  (4,11)  is  insignificant  relative  to  the  others.  Thus: 
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tan6 

tan6„  =  — tan6_  =  tanfi 
E  120  G  m 

p  =  1.78  X  10  "*  ilb(mass)/in® 

It  has  been  assumed  for  simplicity  that  real  parts  of  complex  moduli  and 
loss  tangents  are  frequency  independent. 

For  the  first  mode: 

=  1480  1/sec 

}  Bending  only 

Oi  =  .308  1/sec 
=  1380  1/sec 

}  Timoshenko  beam 

fii  =  4.44  1/sec 

It  is  seen  that  shear  and  rotatory  inertia  have  a  very  small  effect  on  the 
frequency  but  increase  the  attenuation  by  a  factor  of  14.4.  This  example 
shows  that  for  damping  of  viscoelastic  fiber  reinforced  beams  shear  and  rota¬ 
tory  inertia  are  of  major  importance. 

4.2  Forced  torsional  vibrations  of 
laminated  cylinder. 

A  thin  walled  cylinder  which  is  laminated  through  its  thickness  is  built 
in  at  one  edge  and  is  subjected  to  a  sinusoidal  forcing  torque  at  its 
other  edge.  Each  lamina  is  uniaxially  reinforced  and  has  the  same  material 
properties  with  respect  to  its  material  axes.  The  laminate  is  symmetric 
with  following  lamination  scheme: 

reinforcement  in  generator  direction  (axial)  -  volume  fraction  v 
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reinforcement  in  +  8  direction  -  volume  fraction  Vg 

reinforcement  in  -  6  direction  -  volume  fraction  v_g 

Vj  =  v_g  (4.16) 

For  the  purpose  of  analysis  of  torsional  vibrations,  the  only  effective 
laminate  property  needed  is  the  effective  complex  shear  modulus  G12.  It 
follows  from  (3.11)  that: 

Gi2  =  Ci212  =  ^°^Ci212t^/h  +  ^  ®^Cl2i2t^Q/h  + 

(4.17) 

®^Ci212t_Q/h 

where  t^  is  the  sum  of  the  thickness  of  the  0=0  laminae,  t^g  and  t_g  -  the 
sums  of  the  thicknesses  of  the  +0  and  -0  laminae,  respectively. 

Now: 

t  /h  =  V  (4.18) 

0  o 

"  '.e/h  =  Vg 

and  from  the  last  of  (3.16) 


(4.19) 


Introduction  of  (4.18-19)  into  (4.17)  yields 


G*i2  =  ^°^C4.4V„  +  2^®^C«,4V„  (4.20) 

o  0 

By  the  last  of  (3.16)  and  from  (2,4),  (4.20)  assumes  the  form: 
G*i2  =  G^v^  +  [|(n+K-2l^G^)sin^28+2G^cos^20]Vg  ^^21) 


-  25 


It  has  been  mentioned  before  (Section  2)  that  the  imaginary  parts  of 
n,i  and  k  can  be  neglected  with  respect  to  the  imaginary  parts  of  and 
for  the  usual  fiber  reinforced  material.  Consequently,  the  separation  of 
(4.21)  into  real  and  imaginary  parts  assumes  the  form: 

ctz  =  +  [y(n*^+k^-2Jl^+Gj)sin^2e  +  2gJcos220]Vq 

(4.22) 

G*l  =  G^v^  +  (|G^sin22g  +  2G^cos^2e)Vg  (b) 


In  the  usual  fiber  reinforced  materials  generally 


y(n^+k^-2il^+G^)  >  2G 


R 

A 


II  I 

It  follows  that  the  shear  loss  tangent  of  the  laminate 

12 

tanS*  -  (4,23) 

Gi2 

has  a  maximum  for  6=0  and  decreases  monotonically  to  a  minimum  for  6=45®. 

On  the  other  hand  shear  strength  is  smallest  for  6=0  and  increases  monotoni¬ 
cally  to  a  maximum  for  0=45®.  Therefore,  in  design  for  maximum  damping  it 
is  necessary  to  choose  the  smallest  angle  0  which  complies  with  allowable  shear 
stress. 

Let  the  forcing  torque  at  the  edge  xit=£  be  represented  as 


M  =  M  sinwt 
0 
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and  let  the  amplitude  of  angle  of  twist  (j)  at  xi=il  be  written  Amp  (((>) . 
By  standard  theory  of  viscoelastic  vibrations  with  small  loss  tangents 
[2,3] 


cM 

Amp((|))  =  — ^ 


^sin^ (2a) +sinh^ (26) 
cos (2a)  +  cosh  (2B) 


(4.24) 


where 


pj 


a 


ai£ 

c 


B  =  a3)Z;6* 

2c 

p  -  density 

C  -  Torsional  constant  of  section 
J  -  Polar  moment  of  inertia  of  section 
G*2  -  equ.  (4.22a) 

5*  -  equ.  (4,23) 

Numerical  analysis  has  been  carried  out  for  a  laminate  composed  of  boron/ 

epoxy  laminae  with  both  constituent  volume  fractions  equal  to  .5.  Laminae 

fractional  volumes  in  laminate  are 

V.,  =  .6  v.Q  =  .2  with  e  =  22.50“ 
o  Itj 

Analysis  has  been  performed  in  following  stages: 


(a)  Experimental  results  for  epoxy  matrix  complex  shear  modulus  and  loss 
tangent  as  a  function  of  frequency  have  been  described  by  an  empirical 
formula. 

(b)  This  formula  together  with  elastic  properties  of  fibers  have  been  used 
to  compute  effective  complex  moduli  of  the  uniaxially  reinforced  laminae. 
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as  a  function  of  frequency.  For  this  purpose  results  (2.7-9)  and  other 
formulae  given  in  [2,3]  have  been  used. 

(c)  With  the  aid  of  single  laminae  properties  the  real  and  imaginary  parts 
(4.22)  of  the  effective  complex  shear  modulus  Gfj  have  been  computed  as 
function  of  frequency.  It  should  be  noted  that  in  the  present  applica¬ 
tion  1  indicates  generator  direction  of  cylinder,  and  2  the  direction 
normal  to  generator  and  tangent  to  section  contour. 

^  T 

(d)  The  results  for  Gia  and  Giz  have  been  used  to  compute  (4.24)  as  a  func¬ 
tion  of  frequency  to. 

A  plot  of  such  results  is  shown  in  Fig.  3  for  a  cylinder  of  length  5,=100  in. 
and  thin  walled  circular  section.  It  is  seen  that  the  first  resonance  peak  is 
very  significant  and  may  be  regarded  as  an  elastic  resonance.  However,  the 
damping  of  the  viscoelastic  matrix  becomes  more  effective  with  higher  order 
resonances,  the  fourth  one  being  considerably  reduced. 

5.  Conclusion 

It  has  been  shown  that  complex  moduli  of  uniaxially  fiber  reinforced 
materials  and  of  laminates  of  such  materials,  consisting  of  viscoelastic  matrix 
and  elastic  fibers  can  be  computed  in  straight  forward  fashion.  The  results 
can  be  used  for  analysis  of  structural  vibrations  on  the  basis  of  available 
theory. 

Two  structural  examples  have  been  given  to  assess  the  significance  of 
vibration  damping. 

Many  more  other  interesting  applications  can  be  analyzed  by  the  theory 


which  has  been  presented. 
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